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Abstract 

Recently, a global phenomenological fit to high energy data has included charge symmetry break- 
ing terms, leading to limits on the allowed magnitude of such effects. We discuss two possible 
experiments that could search for isospin violation in valence parton distributions. We show that, 
given the magnitude of charge symmetry violation consistent with existing global data, such ex- 
periments might expect to see effects at a level of several percent. Alternatively, such experiments 
could significantly decrease the upper limits on isospin violation in parton distributions. 
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I. INTRODUCTION 



Charge symmetry is a particular form of isospin invariance that involves a rotation of 
180° about the "2" axis in isospin space. For parton distributions, charge symmetry in- 
volves interchanging up and down quarks while simultaneously interchanging protons and 
neutrons. In nuclear physics, charge symmetry is an extremely well respected symmetry, 
generally obeyed at the level of a fraction of a percent Since charge symmetry is 

so well satisfied at lower energies, it is natural to assume that it holds for parton distribu- 
tion functions (PDFs). Furthermore, there is no direct experimental evidence that charge 
symmetry is violated in PDFs. Recently, the experimental upper limits on parton charge 
symmetry violation (CSV) have been improved, by comparing the F\ structure functions 
for deep inelastic scattering induced by muons and neutrinos js, 0, O, 1^. Until now, all 
phenomenological PDFs have assumed the validity of charge symmetry. But a recent global 
fit of PDFs by Martin et al. /T] included for the first time the possibility of charge symmetry 
violating PDFs for both valence and sea quarks. This provides us with parton distribution 
functions that agree with all of the experimental information used to obtain global fits to 
PDFs and which incorporate isospin violation. 

The global fit of Martin et al. assumed a particular form for the charge symmetry violating 
PDFs, for both valence and sea quark CSV. In Sect.^l we will review the form for the CSV 
terms used by MRST, and we will discuss some of the features of their amplitudes. The 
allowed variation in the CSV terms is rather large. Consequently, there may be possibilities 
either to measure the magnitude of parton CSV distributions, or to provide more strict 
experimental limits on these quantities, with dedicated experiments that focus on specific 
observables. We will review two such possibilities in this paper, both of which are sensitive 
to valence quark CSV. 

One promising approach would be to measure Drell-Yan cross sections induced by charged 
pions on an isoscalar target, for which the simplest would be the deuteron. In Sect. IIIH we 
review the possibilities for such measurements. We describe the observable which would be 
most sensitive to valence parton CSV distributions, and we use existing parton distributions, 
plus the CSV terms of MRST, to show potential variations in these observables which 
correspond to the current limits on CSV obtained in the MRST global fit. 

A second possibility is through semi-inclusive deep inelastic scattering for electrons on 
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isoscalar targets. In Sect. IIVI we describe the relevant observable and show the magnitude 
of the effects which one must measure in order to establish charge symmetry violation in 
parton distribution functions, or to provide more strict upper limits on CSV terms. 

We then review the information that could be extracted from these experiments, and the 
prospects for such measurements. 



II. PHENOMENOLOGICAL CHARGE SYMMETRY VIOLATING PDFS 



Because CSV effects are extremely small at nuclear physics energy scales [1|, it is 
natural to assume that parton distribution functions obey charge symmetry. Furthermore, 
there is no direct evidence for violation of parton charge symmetry ^, even though the 
existing direct upper limits on parton charge symmetry violation are at roughly the 5-10% 
level jsl. Martin, Roberts, Stirling and Thorne (MRST) have studied the uncertainties 
in parton distributions arising from a number of factors, including for the first time isospin 
violation. 

Charge symmetry violating PDFs involve the difference between, say, the down quark 
PDF in the proton and the up quark in the neutron; thus we define 

5u,{x) = - , (1) 

with analogous relations for antiquarks. Now, from valence quark normalization, the first 
moment of the valence quark CSV PDFs must vanish, i.e., 
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C dx5d,{x) = f dx - = , (2) 

I' 

with an analogous relation for the first moment of 5u^{x). 

The MRST group chose a specific model for valence quark charge symmetry violating 
PDFs. They constructed a function that automatically satisfied the quark normalization 
condition of Eq. namely: 

5u^{x) = —Sd^(x) = Kf{x) 
f{x) = (1 - (x - .0909) . (3) 

The function f{x) in Eq. Q was chosen so that at both small and large x, f{x) has a form 
similar to the MRST2001 valence quark distributions [^l, and the first moment of f{x) is 




FIG. 1: profile for phenomenological isospin violating parton distributions, for sea quarks 
(top curve) and valence quarks (bottom curve), from the MRST group, Ref. The quantity S 
associated with sea quark isospin violation is defined in Eq. (0J), while the coefficient k is defined 
in Eq. ©. 

zero. The functional form of the valence CSV distributions guaranteed that Su^ and Sd^ 
would have opposite signs at large x, in agreement with theoretical models for parton CSV 
Q,^^. Fixing the form of the CSV parton distribution leaves undetermined only the overall 
coefficient n, which was then varied in a global fit to a wide range of high energy data. 
The value of k, which minimised the in the MRST global fit was k = —0.2. The MRST 
vs. K is shown as the bottom curve in Fig. ^ Clearly has a shallow minimum with 
the 90% confidence level obtained for the range —0.8 < k < +0.65. The top figure in Fig. |21 
plots the valence quark CSV PDFs corresponding to the MRST best fit value, n = —0.2. 
The best-fit phenomenological valence CSV PDFs look extremely similar to the CSV PDFs 
calculated by Rodionov et al. who implemented charge symmetry violation in simple 

quark models; this is shown as the lower figure in Fig. |21 Note, however, that within the 
90% confidence region for the MRST global fit, the valence quark CSV could be up to four 
times as large as that predicted by Rodionov et al. and it could even have the opposite sign. 

The magnitude of the allowed CSV effects obtained by MRST is consistent with un- 
certainties in phenomenological PDFs. The total momentum carried by valence quarks is 
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FIG. 2: Left: the phenomenological valence quark CSV function from Ref. corresponding to 
best fit value k = —0.2 defined in Eq. Solid curve: x6dv{x); dashed curve: xSu^ix). Right: 
theoretical CSV PDFs by Rodionov et ai, Ref. ^]|. Solid line: x6uv{x)] dash-dot line: x6dv{x). 

determined to within about 2%. Inclusion of a valence quark CSV term changes the mo- 
mentum carried by valence quarks in the neutron from those in the proton. The total 
momentum carried by valence quarks is given by the second moment of the distribution, 
e.g., the momentum carried by up valence quarks in the neutron 

"1 



The 90% confidence limit in the valence quark CSV terms corresponds to a variation of 
roughly 2% of the momentum carried by valence quarks in the neutron, or just about the 
known experimental uncertainty in this quantity. 

At the level allowed by MRST, isospin violating PDFs are sufficiently large that, by 
themselves, they could account for the entire anomaly in the Weinberg angle suggested by 



the NuTeV experiment 



12 



M At present, tins .s the only smgle effect that appears 
capable of removing 100% of the NuTeV anomaly [16[; consequently, it is quite important 
that one be able to test the magnitude of parton isospin violation. 

The MRST group also searched for the presence of charge symmetry violation in the sea 
quark sector. Again, they chose a specific form for sea quark CSV, dependent on a single 
parameter. 



M-(a;) = dP{x) [1 + S] 
r{x) = u^ix) [I - 5] 



(4) 
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With the form chosen, the total momentum carried by antiquarks in the neutron and proton 
are approximately equal. 

Perhaps surprisingly, evidence for sea quark CSV in the global fit was substantially 
stronger than that for valence quark CSV. As shown in the top curve in Fig. Q the best 
fit is obtained for S = 0.08, corresponding to an 8% violation of charge symmetry in the 
nucleon sea. The corresponding to this value is substantially better than with no charge 
symmetry violation, primarily because of the improvement in the fit to the NMC fi — D DIS 



data [5|,|6| when -u" is increased. The fit to the E605 Drell-Yan data [17 1 is also substantially 
improved by the sea quark CSV term. 

As explained in the following sections, we have used the MRST CSV parton distribu- 
tions to calculate the differences one could expect in observables for pion-induced Drell-Yan 
processes, and semi-inclusive charged pion production in electron-deuteron deep inelastic 
scattering. 

III. PION-INDUCED DRELL-YAN PROCESSES AND PARTON CHARGE SYM- 
METRY 

In Drell-Yan (DY) processes [l^, two hadrons collide at high energies, and a quark in 
one hadron annihilates an antiquark of the same flavor in the other hadron, producing a 
virtual photon which subsequently radiates a pair of muons with opposite sign. This is 
shown schematically in Fig. El for A^A^ and vrA^ DY processes. 

Since the PDFs for nucleons are rather accurately known, DY processes for pions on 



nucleons provide a sensitive way o: 



extracting parton distributions in the pion. For example. 



NAIO exp.i.e. a. CERN expen^en. a. E..iLab Q bo.. ..died 

Drell-Yan reactions induced by vr^. Such reactions could be studied using the Main Injector 
at FermiLab to produce high energy protons that are scattered from a nuclear target to 
produce charged pions, and subsequently observing Drell-Yan reactions for these pions on 
an isoscalar target, which we will assume is the deuteron. If this program were carried out, 
then from the DY reaction the pion PDFs could be very accurately determined at this value 
of Q^. We will show that by comparing DY cross sections for Tr"*" and tt^ on the deuteron, 
one can extract the CSV violating part of the nucleon valence PDFs - alternatively, through 
such reactions one could produce considerably stronger upper limits on the parton CSV 
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FIG. 3: Schematic picture of the Dreh-Yan process; a quark and antiquark of the same flavor 
annihilate to form a virtual photon that decays into a high-mass muon pair. Left: In NN DY 
processes, a quark in one nucleon annihilates with an antiquark in the other nucleon. Right: in 
TT"*" — p DY, a valence d in the tt^ can interact with a down quark in the proton. 

terms. We use the nucleon CSV distributions recently extracted by MRST to indicate the 
range of CSV that could be probed in such experiments. 

To test nucleon valence quark CSV, it is necessary to measure n — D DY processes at 
kinematics corresponding to large x for both pion and nucleon, e.g., x,x-,t > 0.3. In this 
region, where it is a reasonable first approximation to neglect sea quark effects, the Drell-Yan 
process will predominantly occur when a valence quark in the deuteron is annihilated by a 
valence antiquark in the pion. Then the lowest order (LO) DY cross sections for — D 
and TT^ — D have the approximate form, 

^^^D - l [dnx)+d^ix)]d-\x^) , 

^"-D ^ I [u^{x)+u"ix)]u--{x^) . (5) 

Now, from charge conjugation invariance and the assumption of charge symmetry for 
pion PDFs, we can write the relations, 



TTv{x) 


= ufix) 


= df 


{x) = d:'ix) = ui~ix) 


5 


7rs(x) 


= qfix) 




[x) = qf{x) = qf{x) 


[q = u, d] 


TTs{x) 


= s-'^ix) 




{x). 





(6) 

Inserting the results of Eq. (jH)) into Eq. (0), we see that in the valence-dominated region for 
both pion and deuteron, the tt^ — D DY cross section should be four times the ir'^ — D term. 



Thus, in the quantity 

4a^T^ - a-I^ , (7) 

the valence-valence contributions will cancel. The remaining terms will contain sea-valence 
interference terms and charge symmetry violating terms for both nucleon and pion. Intro- 
ducing the nucleon CSV terms into the PDFs, the DY cross section for tt'^ — D has the 
form, 

9o"^+£,(x^, x) = T[^{xTr) dv{x) + Uv{x) — 6uv{x) + 5 (u{x) + d{x)^ — 6u{x) — A6d{x) 



5 {dv{x) + Uv{x)) + 10 (u{x) + d{x)^ — 6uv{x) — A6dv{x) 



- 26u{x) - 86d{x) 



2tt,{x^) [s{x) + s{x)] , (8) 



in Eq. (jHJ we have also introduced the charge symmetry violating PDFs from Eq. (jT]). 

As proposed by Londergan et al. 21], once the DY cross sections for vr"^ and vr" on 
deuterium are measured, one can extract values for the pion valence and sea quark parton 
distributions. Using those distributions, one can then focus on the region of large x for 
both pion and nucleon. In this region, one can construct ratios of the DY cross sections, for 
example 

R^^ix^, x) = -^"^ (9) 

^n' D " 7r+ D 

As we will see, in order to extract CSV terms in this ratio, the large quantities in the two 
terms in the numerator of Eq. Q will very nearly cancel. This requires that one be able 
to obtain very accurate relative normalization of DY cross sections for charged pions. This 
can be achieved by normalizing the charged pion cross sections to the J/ip peak. 

Inserting the DY cross sections from Eq. (jH)) (and the analogous tt~ —D term), one obtains 

RnEi^n,x) ^ ^^^^ + C{x^)[Ras{x) + R,Ax)]+Rs{x^,x) 



Rcs{x) 
Rsv{x) 
Rsix-Kj x) 



A{6d^{x) - 6uy{x)) 
3q^{x) 

5{u{x) + d{x)) — 6u{x) — 4:6d{x) 
Qvix) 

2Tfs{x^) {s{x) + s{x)) 
t!-v{xt,) q^{x) 
27rs(x^)\ 



CM = 1 , ... 

V -n^ix^) J 

q^{x) = u^{x) + d^{x) (10) 



Eq. (fTUI) is calculated to lowest order in small terms. With the exception of the CSV 
terms, all of the parton distribution functions should be known (the pion PDFs at that 
could be extracted from the same Drell-Yan experiment). Our proposal is that one 
construct the ratio R^q and measure its x dependence for fixed x^^ > 0.3. In this case the 
first term in Eq. (|iup is a constant, and that constant will decrease as increases. From 
the known nucleon and pion valence PDFs, one can predict the value of this ratio assuming 
charge symmetry. Deviations of the ratio from this prediction would be evidence for CSV 
in nucleon parton distributions. The terms Rcs{x) in Eq. ()10j) contains contributions from 
valence quark CSV. Although the CSV term is small, with increasing x it will become a 
progressively larger fraction of the ratio, since all other terms are proportional to sea quark 
distributions that fall rapidly with increasing x. 

We used phenomenological parton distributions for nucleons and pions, that could be 
evolved to the region of interest for Drell-Yan processes. To obtain their CSV parton 
distributions, MRST started with the MRST2001 set of PDFs and varied these using a global 
fit, to obtain the best distributions including charge symmetry violation. The resulting PDFs 
differ slightly from the best-fit MRST2001 distributions 1 221. The value of Aqcd and the 

n 

gluon distributions are essentially identical to those in the MRST2001 set [9|]. The input 
parameters for the nucleon PDFs, allowing for valence quark CSV, are 



xu^{x) = 0.129x°-22^(l 

xdJx) = 0.0163x°-2^Hl - a;)'-'' 



1 + 4.89a;°-^ + 69.86a; 



1 + 123.62;°-^ + 76.04X 



xS{x) = 0.215x"''-''^(l - xy-'^ [1 + 3.34x"-' + 11.80xJ (11) 

MRST define the quark sea from the function S{x) using the relations, 

u{x), d{x), s{x) = 0.2S-0.5A, 0.2S + 0.5A, OAS 
xA{x) = x{d{x) — U{x)) 

= 1.195x^-2^(1 -x)^-^° [l + 14.05a; -45.52a;2] (12) 

The starting scale for these PDFs is Qq = 1 GeV^, and MRST provide interpolation matrices 
that allow one to evolve these PDFs to higher Q^. For the valence quark CSV we used the 
MRST form of Eq. (jH)), and we varied the overall parameter k, between the range —0.8 and 
+0.65 corresponding to the 90% confidence limit obtained by MRST. We have used the same 
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MRST CSV distributions given in Eqs. Q and (0)) at all Q^, in agreement with MRST, who 
did not include any dependence in the CSV terms in their global fits. 

In the MRST global fits, the valence and sea quark CSV were varied separately, and 
slightly different minima were obtained for the quark PDFs in these two cases. Thus there 
is some inconsistency in our using sea quark CSV terms with the quark PDFs appropriate 
for the valence CSV global fit. However, the contribution of sea quark CSV to the ratio 
R^D Eq. Q is extremely small since we restrict our attention to the region where both 
X and are large. In this region the inconsistency in our procedure has only a very small 
effect. 

For the pion PDFs we took the parton distributions of Sutton et al. j^], extracted from 
pion Drell-Yan and prompt photon experiments 20, 24|. The Sutton analysis obtained 
different pion PDFs, depending upon the amount of the pion momentum assumed to be 
carried by the sea. We used those PDFs that correspond to the sea carrying 10% of the pion 
momentum. These PDFs are defined for a starting scale Qq = 4 GeV^ and can be evolved 
upwards using interpolation matrices provided by the Durham group For an actual 

71 — D Drell-Yan experiment, the pion PDFs would be extracted directly from the Drell-Yan 
data at that Q^. 

We have plotted the quantity R^^ defined in Eq. (fTUj) as a function of x, for fixed x-,,, at 
= 25 GeV^, which is a representative value of that could be obtained in fixed-target 
Drell-Yan experiments with charged pions at FermiLab. In fixed-target experiments the 
detector configurations preferentially detect results corresponding to large | 
The top graph in Fig. IH plots R^]^ vs. x for = 0.4, while the bottom graph shows the 
same quantity for = 0.8. 

The solid curve in Fig. |3] corresponds to no CSV contribution; the dashed and dotted 
curves include the limits of the phenomenological CSV contributions at the 90% confidence 
level as extracted by the MRST group. The nucleon PDFs are extracted from global fits, 
and the pion PDFs would have previously been measured in this DY experiment. The 
difference between the sohd curve and the two dashed curves gives the magnitude of CSV 
effects (almost entirely nucleon valence CSV) allowed by the MRST phenomenological fit. 
Note that at reasonable kinematic values, the differences between the CSV contributions are 
substantial. At x = 0.4 the two extremes of the CSV contributions are 30-50% smaller or 
larger than the value corresponding to no CSV contribution, while for x ~ 0.7 the difference 
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FIG. 4: The ratio defined in Eq. ^ clS Oil function of x, for = 25 GeV^. Top graph: 

= 0.4; bottom graph: = 0.8. Sohd curve: no CSV terms, k = 0; dashed curve: k = +0.65; 
dotted curve: k = —0.8. 

between the terms is more like a factor of 2. 

The magnitude of the sohd curve gives the precision necessary in the DY measurements 
in order to extract meaningful information. For x^^ = 0.4, the two DY cross sections are 
predicted to cancel to within about 10%. Therefore the charged pion DY cross sections need 
to be measured to a few percent. For = 0.8 and x ~ 0.7, the two DY cross sections 
cancel almost completely with the pion PDFs used in these calculations. However, the 
CSV contributions at the 90% level would allow values between roughly —0.2 and +0.2, so 
measurements of DY cross sections at the few percent level would allow one to discriminate 
between these limits. 

There is an additional contribution to the Drell-Yan ratio of Eq. arising from the 
possibility of charge symmetry violation in the pion PDFs. The largest contribution would 
arise from differences between the valence d in the and the u in the tc~. The additional 
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contribution to the DY ratio has the form 



(13) 



The contribution from pion CSV is a function only of x^. For fixed it will contribute a 
constant amount to the DY ratio. This was estimated by Londergan et al. 21], by taking 
quark mass effects into account in a Nambu-Jona Lasinio model that was used to calculate 
pion valence quark distributions. The charge symmetry violating pion PDFs were calculated 
at a small value of and then evaluated at higher values using DGLAP evolution. 

The predicted effects from pion CSV were quite small. For Xj, ~ 0.4, the effects were 
almost zero, while for x,r ~ 0.8 the contribution from Eq. (fT!?|l was less than 0.01 {2^. This 
term is added to all of the curves in Fig. |31 and increases slightly the uncertainty in nucleon 
CSV that can be extracted from these ratios. 

IV. SEMI-INCLUSIVE ELECTRON-DEUTERON REACTIONS AND PARTON 
CHARGE SYMMETRY VIOLATION 

Another possibility to measure charge symmetry violating terms in parton distribu- 
tions arises in measurements of semi-inclusive charged pion production from lepton DIS 
on isoscalar nuclear targets. Considering semi-inclusive electroproduction of a hadron h 
from a nucleon A^, e + N h + X , the yield of hadron h in such processes is given by 



In Eq. ()14j] . N^^{x,z) is the yield of hadron h from the semi-inclusive DIS (SIDIS) elec- 
troproduction on nucleon A^, qf{x) is the parton distribution for flavor i in the nucleon. 
Similarly D^{z) is the fragmentation function for a quark of flavor i to fragment into hadron 
h, where z is the fraction of the hadron energy carried by the quark. Semi-inclusive produc- 
tion of a charged hadron from a proton thus has the form 



Ar^'^(x,z) = Ee?gf(x)Df(.) 



(14) 



N'P^{x,z) 



= ^[u^{x)Di{z) + u^{x)D^{z)_ 

+ I [d^{x)D^{z) + dP{x)D^{z) + s^{x)Df{z) + 7s^{x)Df 



(15) 



Charge conjugation invariance requires that 



Dt{z) = DZ{z), D^{z)=Dj{z) . 



(16) 
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The above equations hold for electroproduction of any charged hadrons. In the remaining 
discussion we look specifically at electroproduction of charged pions. If we assume charge 
symmetry for fragmentation functions we have the additional relations 

Df{z) = Dfiz), Df{z) = Dl{z) . (17) 

Thus, under the assumption of charge symmetry, the fragmentation of light quarks to charged 
pions can be written in terms of only two independent fragmentation functions, which are 
defined as "favored" or "unfavored" depending upon whether the quark that produces a 
hadron exists in the valence configurations of that hadron. Thus for up quarks is 
"favored" while D'^ is "unfavored". 

From Eq. (fT^ . the yields depend on "favored" and "unfavored" fragmentation functions, 
as well as the fragmentation of strange quarks to pions. Levelt, Mulders and Schreiber 26| 
derived an expression for measuring the ratio of fragmentation functions. They showed that 

{NP''^ (x, z) - AT"^^ (x, z) + NP""' (x, z) - N""^' (x, z)) _ 95^ A(2) + 1 
(A^?"^+(x,z) - A^'^'^'^(x,z) - A^?"^"(x,2) + iV"'^"(x,z)) ~ 5 A(^) - 1 ^ ' 

In Eq. ()18|1 . the brackets denote integration of both numerator and denominator over x, the 
quantity A (2;) is the unfavored/favored ratio 

and Sa is the Gottfried sum rule j^]. Thus the semi- inclusive electroproduction of charged 
pions m e — D reactions can be used to extract the unfavored/favored ratio of fragmentation 
functions of quarks to pions. 

Londergan et al. j28| showed that these semi-inclusive reactions could also be used to 
investigate the presence of charge symmetry violation in nucleon parton distributions. This 
follows from the fact that, at large x and assuming parton charge symmetry, the favored 
production of charged pions from valence quarks obeys the relation 

iV^"^ (x, z) ^ 4iV^"" (x, z) (20) 

Consequently, they proposed measuring the ratio 

^ 4iV^--(x,z)-iV^-"(x,.) 
^ ' ' - iV^-+ fx, z) - iV^-- (x, z) ^ ' 
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It is convenient to define the quantity 



R^{x,z)^l-^R{x,z) . (22) 



Note that the overall z-dependent factor in Eq. (j22j) is just the factor that is extracted using 
the Levelt et al. ratio defined in Eq. p8|). To lowest order in small quantities it can be 
shown that 

R^{x,z) = Rf{z) + Rcs{x) + Rsv{x) + R^{x,z) 



5A(z) 



^ ^ ^ 1 + A(z) 
oD. . = Df{z) + Dfiz) {s{x) + -s{x)) 

' ^ " ^rW(l + A(z)) qAx) • ^''^ 
The quantities Rcs{x), Rsv{x) and q^ix) are those defined in Eq. (fTn|) . Thus the quantity 
R^{x, z) can be divided into three parts: a term that depends only on z] a term that contains 
the parton CSV contribution, and that depends only on x; and terms that contain contri- 
butions from sea quarks. The term Rf{z) depends on the unfavored/favored fragmentation 
function ratio, which can be accurately measured in these semi-inclusive reactions using the 
method proposed by Levelt et al. (there is a contribution from charge symmetry violation in 



the fragmentation functions; this term was discussed in Ref. j22|. It is expected to be small, 
and it depends only on z). In order to isolate the CSV contribution, one should examine 
the ratio R^{x,z) of Eqs. (PT|) and (f^^ as a function of x at fixed z. Since the CSV term 
is expected to peak at x ~ 0.35 (see Fig. EJ, and the sea quark contributions decrease very 
rapidly at large x, it is possible that at fixed z and sufficiently large x, the CSV terms will 
be substantial, and might even dominate. 

For our calculations we used the fragmentation functions extracted by Kretzer, Leader 
and Christova js^. They used the information on e~^e~ production on charged pions at the 
Z° peak js^, together with the HERMES data on SIDIS charged pion production It 
turns out that the e'^e' data at the peak essentially fixes the combination 



7-)7r+ 



Df + Df + Df 



(24) 



This data, combined with the HERMES SIDIS data, allowed Kretzer et al. to extract the 
individual fragmentation functions, despite the fact that the data need to be evolved 
over a great distance to match up with the HERMES data. They obtained the relations 

Df{z) = 0.689z-i'°=^9(l -2)i-24i 
13 
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FIG. 5: Ratios of fragmentation functions vs. z. Solid line: R^{z); dashed line: z-dependent factor 
in R^{x,z), defined in Eq. I|23|) . Dot-dashed line: the factor C^(z) defined in Eq. (|27jl . Curves 
are calculated for = 2.5 GeV^. 

Df{z) = 0.217z-i-«°^(l - 

Df{z) = 0.164^-1-927^1 _ ^)2.886 (25) 

The fragmentation functions in Eq. ()25|1 are appropriate for an average (Q^) = 2.5 GeV^. 
In this process, the fragmentation functions -DJ^ and D^^ are quite accurately determined, 
while the strange fragmentation function is determined to within about a factor of 2. 

In Fig. we plot the ratios of fragmentation functions that enter into Eq. ()23p. The 
solid curve is R^{z), while the dashed curve is the 2;-dependent factor in the term R^{x, z). 
In Fig. ini we plot the contributions to the quantity R^{x, z) vs. x for fixed z = 0.4. At this 
value, (z) ~ 1.4. The solid curve corresponds to the case with no parton CSV contribution 
(this includes the sea and strange quark contributions). The dashed and dotted curves 
include valence parton CSV contributions corresponding to the MRST phenomenological 
PDFs with K = —0.8 and +0.65, respectively. These values demarcate the 90% confidence 
limits using the MRST valence quark CSV function of Eq. (jSJ. The dot-dashed curve is the 
strange contribution R^{x,z) calculated at z = 0.4. Except for extremely small values of 
X, the sea quark contribution is negligible. Hence the factor 2 uncertainty in the strange 
fragmentation function does not affect this ratio. 

For X > 0.4, the contributions from charge symmetry violating PDFs are substantial, and 
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FIG. 6: Contributions of various terms to the ratio R^{x,z) defined in Eq. (|22|) vs. x at fixed 
z = 0.4. Solid curve: sum of nonstrange and strange sea contributions, Rsv{^) + R^i^:^ = 0.4). 
Dash-dot hne: strange sea contribution, calculated for z = 0.4. Long dash-dot (dotted) curves: 
inclusion of nucleon CSV terms from MRST global fit defined in Eq. H10() . for k = —0.8 and 
K = +0.65, respectively. Curves are calculated for = 2.5 GeV^. 

they rapidly become the dominant contribution at larger x. Thus, at the levels determined 
by the MRST global fit, it would appear that precise measurements of charged pion produc- 
tion in semi-inclusive DIS electroproduction reactions on deuterium have the possibility of 
observing these isospin-violating effects. At the very least they would be able to lower the 
current allowed limits on partonic CSV effects. 

Note that the expected effects shown in Fig. IHl are small. They would require extremely 
precise measurements, in order to be able to distinguish between various predictions for 
valence quark CSV terms. In addition, Eq. shows that the x-dependent terms plotted 
in Fig. inisit atop an x-independent term. For z = 0.4 this term has a value of approximately 
1.4; hence this term is much larger than the CSV term one wants to extract. This difficulty 
could be overcome, provided that one can measure very accurate values for the fragmentation 
functions. In this case one can construct the ratio 




8 



( 



AT^^ {x,z) _ N^"" {x,z) 
l+4A(z) 4 + A(z) 



) 



(26) 



N^^^{x,z) - N^^-{x,z) 
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It is straightforward to show that 



R^{x, z) = C^{z) [Rcs{x) + Rsvix) + R^{x, z) 

C^(z) - ^(^ + ^(")) (27) 
^ - (1 + 4A(^))(4 + A(^)) ^^^^ 



Eq. (jzoj) has the advantage that it ehminates the large 2;-dependent term in Eq. fl23|) . How- 
ever, uncertainties in the favored/unfavored fragmentation ratio now play an important role 
in determining the uncertainty associated with extracting the charge symmetry violating 
PDFs. The term C^{z) in Eq. (P7|) is plotted as the dot-dashed curve in Fig. El at = 2.5 
GeV^. It is normalized so that ~ 1 for moderate values of z. 

Using either of the ratios defined in Eqs. (PHj) or ()26|) . in order to test for CSV terms 
in electron induced SIDIS reactions it is essential that one be able to vary x while keeping 
z constant. Furthermore, the ability to write the ratio as the sum of terms in x and z 
requires that the cross section factorize as in Eq. (fT^ . Thus the data must be taken at 
sufficiently high energies that factorization is valid to within a few percent. SIDIS charged- 
pion experiments have been carried out at HERMES. Another possibility is measurements 
of e + D ^ vr^ -|- X at Jefferson Laboratory. However, 6 GeV is most probably too low an 
energy for factorization to be valid, and even 12 GeV may not suffice. It will be necessary 
to perform checks of the validity of the factorization hypothesis. If an electron-ion collider 
were built, these experiments could be carried out at sufficiently high energies. 

The CSV effects shown in Fig. |B] are slightly underestimated. The relative magnitude of 
CSV and sea quark terms is the same for the Drell-Yan ratio and electroproduction, as can 
be seen by comparison of Eqs. ()10|) and ()23p . The DY cross sections are evaluated at much 
larger values; at these values the sea quark distribution is shifted to significantly smaller 
X. However, we have not evolved the CSV distributions, in agreement with the procedure 
used by MRST in their global fit. If we evolved the CSV distributions in Q^, this would 
increase the CSV contributions relative to the sea in Fig. IHl Finally, our results have been 
derived in lowest order QCD; it is important to check how these results change in NLO. 



V. CONCLUSIONS 



The MRST phenomenological global fit to parton distribution functions provides limits 
on the magnitude of isospin violating PDFs. The CSV effects allowed by the MRST fit 
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are substantially larger than theoretical predictions of charge symmetry violation in parton 
distributions [10|, In this paper, we have analyzed two dedicated experiments that 

might detect parton CSV effects. We used the range of values allowed by MRST to assess 
the magnitude of effects that might be expected in Drell-Yan reactions induced by charged 
pions on the deuteron, and in semi-inclusive production of charged pions in electron-deuteron 
deep inelastic scattering. 

In the Drell-Yan reaction one compares the cross sections for — D and —D reactions 
at reasonably large values of x for both pion and nucleon. Using parton distributions for 
both nucleon and pion, we predict quite large values for this ratio, depending on the sign 
and magnitude of the CSV distributions. 

In the electron SIDIS measurements, once again one compares rates for production of 
7r+ and tt^. In order to extract CSV terms in SIDIS reactions, it is necessary to make very 
precise measurements of both the cross sections and fragmentation functions. It is essential 
that the factorization hypothesis for the cross section be valid to a few percent. Nevertheless, 
given the importance of these quantities and their relevance in experiments like the NuTeV 
neutrino cross sections, it is of great interest now to investigate this issue experimentally. 
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